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^ 1 Abstract 

We have developed a mathematical model of regulation of expression of the Escherichia coli lac operon, 
and have investigated bistability in its steady-state induction behavior in the absence of external glucose. 
Numerical analysis of equations describing regulation by artificial inducers revealed two natural bistability 

> 

' parameters that can be used to control the range of inducer concentrations over which the model exhibits 

CN ■ bistability. By tuning these bistability parameters, we found a family of biophysically reasonable systems 

. that are consistent with an experimentally determined bistable region for induction by thio-methylgalactoside 

o 

(Ozbudak et al. Nature 427:737, 2004). The model predicts that bistability can be abolished when passive 



transport or permease export becomes sufficiently large; the former case is especially relevant to induction 



^ . by isopropyl-/?, D-thiogalactopyranoside. To model regulation by lactose, we developed similar equations 

in which allolactose, a metabolic intermediate in lactose metabolism and a natural inducer of lac, is the 
inducer. For biophysically reasonable parameter values, these equations yield no bistability in response 
to induction by lactose; however, systems with an unphysically small permease-dependent export effect 
can exhibit small amounts of bistability for limited ranges of parameter values. These results cast doubt 
on the relevance of bistability in the lac operon within the natural context of E. coli, and help shed light 
on the controversy among existing theoretical studies that address this issue. The results also suggest an 
experimental approach to address the relevance of bistability in the lac operon within the natural context of 
E. coli. 
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INTRODUCTION 



In 1957, Novick and Weiner discovered that Escherichia coli can exhibit discontinuous switch- 
ing in expression of the lac operon, with some cells expressing a large amount of /3-galactosidase 
(/3-gal), other cells expressing a small amount, and an insignificant number of cells expressing an 
intermediate amount [1]. Recently, this effect was further characterized using single-cell assays of 
fluorescence levels in a population of E. coli cells carrying a lac::gfp reporter [2]. The population 
exhibited a bimodal distribution, with induced cells having over 100 times the fluorescence level 
of uninduced cells. These observations have been attributed to the existence of two steady states, 
i.e., bistability, in the induction of lac in E. coli. 

Recent modeling studies have emphasized the importance of determining whether bistability in 

n n n n n 

expression of lac is relevant within a natural context O, |4j, [5|, [6|, |7J]. This question remains open 
because experimental studies have focused on the response of lac expression to artificial inducers, 
such as thio-methylgalactoside (TMG) and isopropyl-/?, D-thiogalactopyranoside (IPTG), rather 
than the natural inducer, allolactose. This difference is critical because artificial inducers (also 
known as gratuitous inducers) are not metabolized by the induced enzyme, whereas the natural 
inducer is a metabolic intermediate in lactose degradation, which is catalyzed by the induced 
enzyme. 

Savageau [3] found important differences between induction by IPTG vs. lactose in his theoret- 
ical treatment of bistability in the lac operon. In Savageau's model, because production and decay 
of allolactose are both proportional to the /3-gal concentration, bistability is forbidden. Expression 
of lac in response to lactose was therefore predicted not to exhibit bistability. This prediction 
agreed with the absence of steady-state bistability in an experimental study of populations of E. 
coli cells exposed to lactose, described in the Supplementary Material of Ref. D2D — in that study, 
only transient bimodal distributions of green fluorescence levels among cells were observed at 
some glucose concentrations. It was later noted that models with operon-independent decay of 
lactose (e.g., due to dilution by cell growth) could exhibit bistability tn\- Several studies using 



such models found either a bistable or graded response to lactose, depending on parameter values 
or external glucose levels [5,^ 



L il 7, 8, 9J], and in agreement with the model of Savageau, a model of 



van Hoek & Hogeweg QVQ was explicitly shown to exhibit no bistability in the absence of operon- 
independent decay of allolactose. However, these studies disagree in their assessment of whether 
bistability is present |5l Is |9J] or absent [7|, [8Q in expression of lac among E. coli cells in a natural 
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context. 

In addition to predicting whether lac induction exhibits bistability, some studies have addressed 
the question of whether bistability might enhance or hinder the performance of E. coli cells. Both 
Savageau [4] and van Hoek & Hogeweg [8] found that bistability increases the time required to 
respond to sudden increases in environmental lactose, which can be a disadvantage in competition 
for nutrients. These results argue against the natural relevance of bistability in lac expression. 

Another important question that has not yet been addressed is whether the experimental obser- 
vations of bistability in Ref. [|2|] are consistent with independent biophysical data that characterize 
processes relevant to regulation of lac expression. Although phenomenological models were de- 
veloped to reproduce the steady-state behavior [|2|] and the experimentally characterized dynamics 



of switching between stable steady states [10], these models were not constrained by independent 
biophysical data. For example, it is unclear whether the phenomenological models are consistent 
with independently measured permease transport kinetics. On the other hand, studies of bistability 
using more detailed, biophysical models of lac induction were either only partially constrained [|7|] 



or did not consider the response to artificial inducers [|5L 

Here we analyze bistability in an ordinary differential equation (ODE) model of lac induction. 
We use ODEs because we restrict our analysis to steady-state behaviors, and because the protein 
concentrations in fully induced cells are O(10 4 ) per cell (see Parameter Values section) and have 
negligible fluctuations. Similar equations describe induction by lactose or artificial inducers. We 
first use the model to gain insight into key determinants of bistability of lac expression in response 
to artificial inducers, and to understand how characteristics of bistability are controlled by model 
parameters. We then use the resulting insight to tune the parameters of the model to match the 
bistable behavior observed by Ozbudak et al. (20, and to predict mechanisms by which bistability 
might be abolished. Finally, like previous modeling studies, we use the model to address the 
question of whether lac expression might be bistable in a natural context, contributing to resolution 
of what is now a long-standing controversy. 

MODEL 

In our model of lac induction (Fig.Hk), the following set of coupled ordinary differential equa- 
tions relate the internal lactose concentration (/), allactose concentration (a), and /3-galactosidase 
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FIG. 1: Circuitry for models of lac induction, a) Model for induction by lactose (Eqs. (Q~|)), including 
the following processes: (1) proportional production of permease (LacY) and /3-gal (LacZ); (2) permease- 
mediated transport of lactose; (3) dilution of intracellular species by cell growth; (4) /3-gal catalyzed degra- 
dation of lactose, producing both the metabolic intermediate allolactose, and the ultimate products of degra- 
dation, glucose and galactose; and (5) /3-gal catalyzed degradation of allolactose, producing glucose and 
galactose, b) Model for induction by artificial inducers (Eqs. Q), including: (1) proportional production of 
permease (LacY) and /3-gal (LacZ); (2) permease-mediated transport of inducer; (3) dilution of intracellular 
species by cell growth and (6) passive transport of inducer. 



concentration (z) to the external lactose concentration (/*) 



(I* ~ g) Pzl 

OLZ — — — — Tt (la) 

Ki + l + l* (l + a/K m>a )K m>l + l 
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(1 + a/ K m>a ) K m> i + 1 (1 + 1/ K m j) K m>a + a 
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In Eqs. ([T]), a and (pa are the rate constants for permease-dependent lactose import and export, 
respectively, i<Q is the Michaelis constant for permease-dependent lactose transport (assumed to 
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be the same for import and export), (3 and K m i are the rate constant and Michaelis constant for 
lactose degradation, v is the branching fraction of lactose degradation to allolactose, 5 and K m>a 
are the rate constant and Michaelis constant for allolactose degradation, 7 is the rate of dilution 
due to cell growth, cy and are the basal and inducible rates of /3-galactosidase production, K z is 
the allactose concentration at half-maximal induction of /3-galactosidase production, and n is the 
Hill number for lactose induction of /3-galactosidase production. 

The metabolic fluxes in Eqs. (OQ) include the effects of competition between allolactose and 
lactose for access to /3-galactosidase (/3 and 5 terms). Because shuttling of galactosides across 
membranes occurs through a single permease channel [11], we also consider the influence of 
competition between external and internal lactose for access to permease (a and <p terms); however, 
as a simplification, we do not consider transitions among distinct internal states of the permease 

y. 

To focus on the operating conditions of the system that are most relevant to lactose utilization by 
E. coli, we only consider regulation in the absence of glucose. This focus is appropriate because, 
in the presence of glucose, lac is not essential for growth, and induced /3-galactosidase levels are 

low nidi. 

Similarly, the model of artificial induction of lac (Fig.[Tb) is given by 



«o (t — £j + ctz— — - i - — 7/ and 



07 + 
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Ki + l + l* 
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In Eqs. ©, variables and parameters have the same meaning as in Eqs. (0Q), except I and I* 
correspond to the level of internal and external artificial inducer (e.g., IPTG or TMG), respectively, 
and a is the rate constant for leakage across the membrane. 

In Eqs. (0Q) and Eqs. ©, protein expression is lumped with gene expression, and the dependence 
of promoter activity on the level of signal (IPTG, TMG, or allolactose) is mode 
Hill function, which is significantly simpler than other models J5, 6, 7, 8l.l9l. fi3i 
hand, Eqs. © considers the effects of competition among substrates in permease transport and 



ed using a simple 
On the other 



strates in permease transport ai 
JflflflflQQ. Compared 



to 



metabolic processes, unlike other models of lac induction [|2j,|3 , 
the model of Savageau Q3|, |4£| , Eqs. (OQ) considers operon-independent decay of allolactose, without 



L3, 



41lZD, as discussed above. Overa 



which bistability in response to lactose is impossible 
and Eqs. © are less detailed than the lac induction models used in Refs. [|13|l 



. Uveraii, hfls 
, H A, fl, 



1, Eqs. © 



and [9], and are more detailed than those used in Refs. |3], fl, Ill5ll . and fl2|], and they therefore 
constitute intermediate complexity equations describing lac induction. Compared to the simpler 
models, the intermediate level of detail provides increased contact between model parameters and 
biophysical measurements, and compared to more detailed models, it facilitates analysis of the 
equations and interpretation of the results. 



PARAMETER VALUES 



We used the parameter values and ranges listed in Table H to analyze bistability in Eqs. (Q~|) and 
Eqs. 0. The values in the table were obtained as follows: 



• 7. We assume the generation time under the conditions in Ref. 120 is 30-60 min. We note, 
however, that this time might be very different for E. coli growing under stress in the gut; this 
represents a source of uncertainty concerning the biological relevance of our predictions. 

• a . We assume that a = except for the case of IPTG, where we explore a range consistent 
with that considered in Ref. [8|]. 

• a. An approximate range of 1-100 s -1 for sugar transport turnover numbers was obtained 
from the review by Wright et al. |l6|] . The range is broader than measured values |r7|] because 
measurements were made at 25°C rather than at the physiological temperature of 37°C in 
the host environment of the gut that we are focusing on here, and at which measurements in 
Ref. [2] were performed. The nominal value of 1000 min -1 was estimated from Ref. [| 1 VFI 
assuming the production rate of permease is the same as that of functional /3-gal. Because 
permease is a monomer while /3-gal is a tetramer, this assumption entails a four-fold smaller 
production rate for permease. This seems possible, as (1) galactoside acetyltransferase 
(GATase) monomer synthesis is eight-fold smaller than /3-gal monomer synthesis; (2) due to 
incomplete operon transcription and the order of genes in the operon (lacZYA), the amount 
of mRNA transcribed from the GATase gene (lacA) and permease gene (lacY) is smaller 
than that from the /3-gal gene (lacZ); (3) there is some evidence that permease is made in 



smaller amounts than /3-gal [18]. 



• 0. We assume no export flux through permease in the artificial induction model, and then 
examine the consequences of introducing such a flux on bistability. Guided by Ref. Ill in . for 
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the lactose model, we use a nominal efflux rate constant {(pa) of half the value of the influx 
rate constant a, and allow the value to decrease in the search for bistable conditions. 



Ki. For simplicity, we assume the same Michaelis constant for permease import and export 



a nominal value of 0.5 mM was obtained from Ref. [17]. The range was applied as per a, 
and encompasses measured values [17, Jjl 2ol . 

• (3. A total lactose turnover number for /3-galactosidase of 2.85 x 10 4 min 1 is estimated 
from a measured value of V max = 61.3 fxmol min" 1 mg" 1 in Ref. 112 ill . This estimate is 
an order of magnitude greater than the value 3.6 x 10 3 min" 1 given in Ref. Il22ll . but the 
two estimates agree closely when one considers that /3-gal converts about half of its lactose 
substrate to glucose and galactose, rather than allolactose, and that the enzyme is composed 



of four monomeric catalytic subunits. The estimate given in Ref. 112211 is appropriate for 
total turnover of lactose on a per monomer basis. Like for a, because measurements were 
performed at 30°C, we consider a range of values ten times lower to ten times higher than 
the nominal value. 



• K m i. The nominal value was obtained directly from Ref. [23]. As for (3, because of 
temperature considerations, we use a range from ten times lower to ten times higher than the 
nominal value. 

• v. The value v = 0.468 was calculated from the total rate of /3-gal degradation of lactose and 
the partial flux from lactose to allolactose reported in Ref. 112 ill . We take it to be a constant 
because the ratio of reaction products was found to be insensitive to temperature changes 
between 30°C and 0°C. 

• 5. An allolactose turnover number for /3-gal of 2.3 x 10 4 min -1 is estimated from a measured 
value of Vmax = 49.6 /x mol min -1 mg" 1 in Ref. |23|]. As for f3, because of temperature 
considerations, we use a range from ten times lower to ten times higher than the nominal 
value. 



K m>a . The nominal value was obtained directly from Ref. II23I1 . As for {3, because of 
temperature considerations, we use a range from ten times lower to ten times higher than the 
nominal value. 
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• e. Using a production rate of 5 /3-gal tetramers per cell per second for a 48 min generation 
time [24], 14,400 molecules are produced during a generation at full induction-this is the 



number of molecules in the cell after doubling (supporting our choice of a noiseless model). 
Assuming a 1 yum 3 mean cell volume ^25\ and linear volume increase in time the 
volume after doubling is approximately 0.7 /xm 3 , leading to a concentration of 34,286 nM. 

c. This value is derived from e, assuming a 1000-fold increase in /3-galactosidase levels 
upon induction Il27ll . 



• K z and n. These values are estimated from IPTG induction data in permease knockout cells 
both from Ref. Il28ll . Fig. 15 and from data compiled in Ref. Il29n . Figs. 1 and 2. The 
nominal value n = 2 was estimated from the slopes of the curves in the figures, and K z 
was determined by estimating from the figures the concentration of IPTG at half-maximal 
induction. The nominal value of 10 5 nM was estimated from data compiled in Ref. 12911 . To 



determine the range, an approximate lower value of 10 4 nM was obtained from Ref. [|28 
and we allowed for an upper value of 10 6 nM to account for potential differences between 
induction by IPTG and TMG or lactose. 



RESULTS 

We first used Eqs. © to determine how parameter values control bistability in the steady-state 
response of lac expression to artificial inducers. To detect and characterize bistability for a given 
set of parameter values, we solved for z(l) and I* (I) as rational functions of /. Bistability in lac 
expression exists when the line describing steady-state levels of z vs. I* adopts a characteristic "S" 
shape, as shown in Fig. [2] Within the bistable range of /*, the highest and lowest levels of z are 
stable steady-state solutions and the intermediate level of z is an unstable steady-state solution of 
Eqs. (OQ). The bistable range is defined by the lower (7* = L) and upper (I* = U) turning points, as 
illustrated in Fig. [2] An analogous signature of bistability can be seen in examining steady-state 
levels of / vs. I* (not shown). For a model with given parameter values, L and U can be located by 
finding the roots of either dl* jdz or dl* jdl using an eigenvalue solver. 

We analyzed Eqs. © for systems with sets of parameter values drawn from the ranges in TablelH 
taking a = 0, = 0, and n = 2. Sets of 100 values each for Ki and K z were obtained using 
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TABLE I: Parameter values. 



Par am 


Description 


Nominal 


Range 


7 


growth rate 




0.0116 min" 1 - 0.0231 min" 1 


a 


passive transport 
rate constant 





- 1.35 min -1 


a 


permease import 
turnover number 


600 min -1 


6 x 10 1 min" 1 - 
6 x 10 3 min" 1 


<P 


ratio of permease 
export to import 
turnover numbers 


(artificial inducers) 
or 0.5 (lactose) 


0-0.5 




permease Michaelis 
constant 


5 x 10 5 nM 


5 x 10 4 nM - 
5 x 10 6 nM 




/3-gal lactose 
turnover number 


2.85 x lO^iin- 1 


2.85 x 10 3 min" 1 - 
2.85 x 10 5 min" 1 


V 


lactose^ allolactose 
/3-gal branching 
fraction 


0.468 




Km 1 


/5-gal lactose 
Michaelis constant 


2.53 mM 


0.253 mM-25.3mM 


5 


/3-gal allolactose 
turnover number 


2.30 x lO^in- 1 


2.30 x 10 3 min" 1 - 
2.30 x 10 5 min" 1 


K 


/3-gal allolactose 
Michaelis constant 


1.2 mM 


0.12 mM -12.0 mM 


e 


fully induced 
/3-gal level 


34285 nM 




c 


basal /3-gal level 


34.3 nM 




K z 


signal level at 
half-maximal 
lac induction 


10 5 nM 


10 4 nM - 10 b nM 


n 


Hill number for 
signal-dependent 
lac induction 


2 





logarithmically even sampling over their allowed ranges. Because the steady-state solutions of 
Eqs. © only depend on a and 7 through the ratio a/7, rather than sampling a and 7 individually, 
we obtained 100 values of a/7 using logarithmically even sampling between the upper and lower 
bound computed from Table HI This sampling scheme yielded 100 x 100 x 100 = 10 6 systems 
with different values of (a/7, Ki, K z ). 

We found that all 10 6 systems exhibited some degree of bistability in response to induction by 
artificial inducers. The dependence of the range of bistability on model parameters was further 
analyzed using two measures that we introduce here: the ratio U/L, and the product U L. We used 



9 




10 

1* (uM) 



FIG. 2: An example of a system from Eqs. ((U) with the upper (U) and lower (L) turning points consistent 
with the results in |2|]. The parameter values are 7 = .0231 min -1 , a = 60 min -1 , K z = 123, 285 11M and 
Ki = 1,077, 217 nM. 



these measures to estimate the percentage of systems for which bistability might be observable in an 
experiment like that in Ref . [2] . By inspecting the measurement errors in Ref. [J, we estimate that 
systems with U/L > 1.1 and UL > 0.01 /iM 2 exhibit bistability that is favorable for experimental 
observation (i.e., difficult to detect), and that systems with either U/L < 1.1 or UL < 0.01 /iM 2 
exhibit bistability that is unfavorable for experimental observation. Among systems with parameter 
values sampled as described above, by these criteria, we predict that experimental observation of 
bistability would be favorable for 65% of systems, and unfavorable for 35% of systems. 

Increasing either a or above zero tends to reduce or abolish bistability in artificially induced 
systems. As a is increased (Fig. [3]), first U begins shifting to lower values of I*, then L begins 
shifting to higher values of I*, leading to an asymptotic behavior in which bistability is abolished. 
Like changes in a , as <\> is increased (Fig.H]), L shifts to higher values of/*; however, by contrast, 
U does not initially show a significant change. As </> is increased further, the entire induction curve 
begins to shift to higher levels of I*. 

* n 

To compare Eqs. d2j) to the data in Ref. B2fl, we first selected a subset of systems for which 
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FIG. 3: Effects of variations in the a>o > parameter on an artificially induced system with = min 1 
and ao = 10~ fc min -1 , k = 0, . . . , 4. The other parameters are given by n = 2, 7 = .0231 min -1 , e = 
34286 nM, c = 34.3 nM, K { = 5 x 10 6 nM, = 10 4 nM and a = 60 min -1 . 

the bistable region is in the same neighborhood as that in Ref. from 3 /iM to 30 /iM TMG. 
Considering this range, out of the 10 6 systems sampled, we selected 187,108 systems for which 
L > 1 /iM and U < 100 /iM for further analysis. Interestingly, we found that all of these systems 
collapse to a single curve when displayed in the space of \og 10 (U/L) vs. \og(Ki/ K z ) (Fig. [5]), 
indicating that U/L can be precisely tuned using the parameter X = Ki/K z . As shown in Fig.[H 
the dependence was accurately modeled using the equation 



log 10 (C//L) 



{Ki/K g 



.9.3 



93 



— > 0. 

10 - 



(3) 



Next, we found that, at a given value of X = Ki/K z , without changing the value of U/L, UL 
could be tuned precisely using the parameter Y = KiK z ^ /a. As shown in Fig. |6j this dependence 
was accurately modeled using the equation 



\og 10 (UL) = C {X) + C X {X) log 10 (F). (4) 
Figure [7] shows the X-dependence of the parameters C (X) and C\{X), obtained numerically 
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FIG. 4: Effects of variations in the > parameter on an artificially induced system with ao = 10 -4 min -1 
and <fi = and 10 _fc min -1 , k = 1, . . . , 4. All of the scales are in fiM. The other parameter values are as in 
Figure [3] 



using systems with similar values of X. For the range of systems considered here, we found that 
Cq{X) could be fit using a third order polynomial in log 10 (X), and that C\(X) could be taken as 
a constant. 

The above phenomenological results provide a prescription for tuning the range of bistability 
exhibited by an artificially induced system. First, the value of U/L can be specified by choosing a 
value of the parameter X = Ki/K z using Eq. ©. Then, using this value of X, the value of UL can 
be specified by choosing a value of the parameter Y = K z Kia/^ using Eq. © and the empirically 
determined Co (X) and C\ [X) (Fig. [7]). We used this prescription to obtain a family of systems that 
are consistent with the parameter values in Table U and that exhibit a range of bistability consistent 
with that observed in Ref. [2|], with log 10 (C//L) w .86 and log 10 (C/L) w 1.92. An example of the 
steady-state behavior of one such system is illustrated in Figure[2] 

We used similar methods to analyze Eqs. ([I]) which describe induction by lactose. No bistability 
was present in the system with nominal parameter values from Table U (Fig. [8]) with = 0.5, which 
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i z 



FIG. 5: Modeling the width as a sigmoid function of Ki/K z . Only a sample of data points are shown. 

is consistent with the theory of Savageau [3] and the Supplementary Material of Ref. 020. However, 
guided by the results for artificial inducers in Fig. HI we examined systems with = 0. Although 
the system with nominal parameter values and = did not exhibit bistability, other systems that 
have parameter values consistent with the ranges in TableUdid exhibit bistability. We then located 
the system that exhibits the largest values of U/L and UL; for this case, a, (3, 5 and K z assume 
their lowest values in Table|I]while 7, K mU K mA and K { assume their highest values (Fig.©. The 
curve in Fig. |9] illustrating bistability characteristics for this system closely resembles a similar 
curve shown in van Hoek & Hogeweg QVU , Fig. 2B. Thus, although our model is less detailed than 
theirs, it can exhibit comparable steady-state behavior. 

To estimate the distribution of systems exhibiting the different qualitative behaviors, as for the 
case of artificial inducers, we analyzed 10 4 systems with randomly sampled parameter values, 
all with = 0. We predict 99.87% of these systems to exhibit no bistability, 0.05% to exhibit 
bistability favorable for observation (U/L > 1.1 and UL > 0.01 /iM 2 ), and 0.08% to exhibit 
bistability that is unfavorable for observation (U/L < 1.1 or UL < 0.01 /iM 2 ). However, as 
observed for Eqs. ©, increasing to even a small fraction of its nominal value rapidly abolishes 
bistability for all combinations of other parameter values in Eqs. © (Fig.fTOl). 
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FIG. 6: Using the Co and C\ in Figure |7] results in a linear relation between the center log 10 (UL) and 
log w (KiK z -f/a). 

CONCLUSIONS 



For the equations describing induction by artificial inducers, we found that the range of external 
inducer concentrations over which systems exhibit bistability is precisely controllable by two 
rational combinations of model parameters. By adjusting these parameters, we were able to 
demonstrate agreement with the bistable range for TMG induction from Ref. (2D. However, in 
achieving this agreement, we assumed that permease-dependent efflux of artificial inducers is 
negligible (0 = 0). We have not found independent biophysical data to constrain this parameter 
for artificial inducers, and therefore predict that it has a value much less than the value of roughly 
0.5 that has been measured for lactose. 

To achieve agreement with the bistable range of roughly 3 fiM to 30 /iM in Ref. |2Q, c and e 
in Eqs. © were tuned to exhibit a 1000-fold induction of protein expression. While this value 
is reasonable based on previous studies, it does disagree with the roughly 100-fold induction of 
GFP expression reported in Ref. [2]. We did analyze systems with alternative values of c and 
e that yield 100-fold induction; however, none of them exhibited bistable ranges that agree with 
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FIG. 7: Dependence of the coefficients in the regression model of Eq. (01) on X = Ki/Kz. 



the range reported in Ref. D2D. Further studies will be required to understand why Eqs. (0 does 
not simultaneously agree with both the bistable range and maximal induction of the experimental 
lac :: gfp reporter system. In addition to model refinement, it would be fruitful to seek systematic 
differences between expression from chromosomal lac and the plasmid-based lac :: gfp reporter 
system used in Ref. D20. 



The lack of bistability observed for induction by lactose agrees with modeling studies concluding 
that bistabity in lac expression is irrelevant to E. coli in a natural context |3, i 3, 8]. Thus, although 
bistable behavior in lac is now well-documented flflfl, because it has only been experimentally 
observed using artificial inducers, its relevance within the natural context of E. coli is doubtful. 
Indeed, it is surprising that the lac operon has been considered to be a paradigm of bistability in gene 
regulation, considering the gaps in understanding that remain after so many careful experimental 
and theoretical studies. 

The present results predict that bistable behavior can be promoted by (1) hindering the kinetics 
of permease transport (a, Kt) and /3-gal catalysis (J3, S, K m j, K m ^ a )\ (2) lowering the required level 
of allolactose for half-maximal lac expression (K z ); and (3) accelerating cell growth (7). These 
predictions suggest genetic targets for engineering E. coli strains that exhibit a clear signature of 
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FIG. 8: Allolactose system with all of the parameters given by their nominal values in Table U 

bistability. Experiments to compare the behavior of such strains with wild-type cells would help 
to clarify whether bistability in lac expression is relevant in a natural context. 
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FIG. 9: Bistability in the = lactose-induced system with a, (5, 5 and K z at their lowest values in TableJ] 
and 7, K m \, K m ^ a and K t at their highest values. This is the system that exhibits the largest values of U/L 
and U L within the allowed ranges of parameter values. 
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FIG. 10: This is the same system as in Figure [9] with <fr = and 10 min 1 , k = 1, . . . , 4. 
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